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Abstract 



We study a model in which p independent Ising spins are coupled to 2d 
quantum gravity (in the form of dynamical planar </> 3 graphs) . Consider- 
ation is given to the p^oo limit in which the partition function becomes 
dominated by certain graphs; we identify most of these graphs. A trun- 
cated model is solved exactly providing information about the behaviour 
■^j- ■ of the full model in the limit f3 — ► 0. Finally, we derive a bound for the 

' critical value of the coupling constant, /3 C and examine the magnetization 

transition in the limit p — > 0. 

o 

On '. 1 Introduction 



Models in which conformal matter of central charge c is coupled to 2d quantum gravity have attracted 
considerable interest. Much progress has been made for the case c < 1, using both continuum and discrete 
approaches. The former has yielded the KPZ formulae jlj |J, which give as functions of the central 
' charge the modification of critical exponents due to gravity. However, these formulae only apply for 

c < 1; for larger values of c they predict that the string susceptibility j s t r is complex, which does not 
make much sense. The discrete methods involve using dynamical triangulations coupled to matter fields. 
Matrix model techniques have proven to be extremely valuable for studying c < 1 models pi J6| and 
perturbative expansions allow one to investigate numerically the behaviour of j str for c > 1 The 
model in which a single Ising spin is attached to the face of each triangle has a central charge of one 
half and has been solved exactly jl0|, [□], |l2], |l3|, , yielding results for the critical exponents that agree 
with those of the KPZ formulae. Recently some remarkable analytical progress has been made for the 
c-^oo limit by considering the low temperature expansion jl5|, [l6| [H]]. There is evidence that in this limit 
tree-like or branched polymer graphs dominate the behaviour of the model and this is supported by our 
results. Despite much effort no one has managed to solve analytically any other models for c > 1; however 
such models are well-defined and there is no obvious pathology at c — 1. Many Monte Carlo simulations 
have been performed jl^, [l9| |2^, ^ 22, ||^| in an attempt to investigate the behaviour near the c = 1 



barrier and to discover whether the breakdown of KPZ theory is due to some change in the geometry, 
such as a transition to a branched polymer phase. So far these simulations have failed to produce any 
convincing evidence of such a phase transition at c = 1. 

In this paper we study the properties of a model for which p independent Ising spins are attached 
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to the face of each triangle, giving a central charge of c = p/2. In section 2 we define precisely three 
slightly different models and in section 3 we examine the limit of large p, identifying most of the dominant 
graphs. A version of the model in which the free energy is truncated is used in section 4 to study the 
transition to behaviour typical of large p, in the limit of small j3. In section || we study the properties of the 
magnetization transition, deriving a bound on the critical value of the coupling constant (/3 C > 0.549) and 
looking at the nature of the transition in the limit p — > 0. In section ^ we conclude by discussing possible 
forms of the phase diagram and by relating our work to the results from various computer simulations 
and to analytical work carried out by other authors. 



2 Definition of the model 

The model is one of p independent Ising spins on each vertex of a random 4> 3 graph (which is the dual 
graph of a triangulated random surface). For a fixed N- vertex 4> 3 graph, G, with a single spin on each 
vertex, the Ising partition function is 

{S} V <ij> J 

with 

Z = 2 N (coshf3)^, (2) 

where Si is the spin on the z-th vertex of the graph, the sum over < ij > is a sum over nearest neighbours 
(by which we mean vertices that are connected by a line in the graph) and (3 is the coupling constant. 
The factor of Zq~ x is introduced in order to simplify the formulae later on. We then sum over a set of 4> 3 
graphs, with N vertices and a fixed genus, g, to get the partition function, 

Z N (p) = J2-(.Z G ) p . (3) 

S c 
G G 

Each graph is weighted by the symmetry factor s G _1 , which is the reciprocal of the order of the symmetry 
group for that graph. The symmetry factors are inserted because they occur in the matrix model solution 
of the p — and p = 1 cases. However, we could equally well take s G = 1 in the above definition, giving us 
a slightly different model, which is nonetheless expected to have identical properties in the thermodynamic 
limit N — + oo. Consequently we will often ignore the symmetry factors especially in the large N limit 
when we would expect the graphs not to be very symmetric and hence to have s G w 1 anyway. 

In this paper we work with planar diagrams (ie g = 0) and consider three different versions of this 
model, which differ with respect to the sets of graphs that are used. In model I, G runs over all the planar 
connected <fi 3 graphs with N vertices. Model II is the same as model I, except that tadpoles are excluded 
so that the graphs are one-particle irreducible. For model III, tadpoles and self-energy terms are excluded 
giving two-particle irreducible graphs. 

In the thermodynamic limit the partition function has the asymptotic form, 

Z N (p) = e »(v,m N y* t r-3 (a Q + ^ + ...y (4) 
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The free energy (i(p, (3) is denned as 

Mp,j8)= Hm ilog(^(p)) (5) 

iv — >oo iv 

and discontinuities in its derivatives indicate the presence of a phase transition. Similarly, fiQ for a given 
graph, G, is defined by, 

M/J)= Hm ibg(Z G ). (6) 

N—>oo iv 

The string exponent 7 str is believed to be universal and to depend only upon certain general characteristics 
of the model. 

2.1 Matrix model results (p = and p = 1) 

The casep = 0, where there are no Ising spins, just corresponds to enumerating <fi 3 graphs. The number of 
graphs Q^(N) in model I can be calculated, for example by matrix model methods Q (see also p5l p6| ) 
with the result, 



y [ ' 2(f + 2)!T(f + 1) ^- • W 

For models II and III, the number of TV- vertex graphs are denoted by Q^ 2 \N) and Q^\N) respectively 
and can also be calculated giving, 

S »(iV) = |^ „~, (9) 

It should be noted that all these models yield the asymptotic form given in (4) with 7 str = — |. 

The p = 1 case has been solved analytically |h], |ll| and has a third order phase transition from a 
disordered to a magnetized state. For model I, the critical value of the coupling constant is given by Jl3[ , 

/3 C = -\ log (± (2V7- lfj « 0.9196, (10) 

and for model III @, 

In both cases, the critical exponents are a = — 1 , j3 = |, 7 = 2, S = 5, ^ = ^- and 77 = 2 — where 
dir is some unknown dimension depending on the geometry of the graphs. Also, 7 str — everywhere, 
except at the critical point, where j* tr = — |. 

These results should be compared with those for a fixed regular lattice, such as the hexagonal lattice, 
for which, 

I3 C = tanlT 1 (Jj^J ~ 0- 658 (!2) 
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and there is a second order magnetization transition with critical exponents: a = 0, /3 = g, 7 = j, S = 15, 
v = \, r\ = \. Introducing the sum over triangulations changes the universality class, but the precise 
nature of the sum is not important. 



3 Dominant graphs in the limit p — > oo 

3.1 Concavity of /j,(p,(3) 

In this section we will ignore the symmetry factors, however it is easy to show that all the results follow 
if we include them. Using the Cauchy-Schwartz inequality on (3) it follows that, 

and hence that (i(p, /?) is concave with respect to p, 

<^Mp,/3)+m(<7,/3)]. (14) 

(A slight modification of the standard argument shows that /j, is also concave with respect to j3). It is 
trivial that for p > 1, 




z N ( P ) = ^(z G r< ( y^z i; j . n.-,) 

so that, 

m(p,/?)<Mi,/?), p>l. (16) 

We also have that, 



N^oo Zm * — ' 
G 

and hence that /i(p, /3) is a concave monotonic increasing function with respect to p. However, fx is bounded 
by a linear function and hence it must be asymptotically linear in p as p — > 00. Since the number of graphs 
with given Zq does not depend upon p it must be the case that, as p — > 00, Z^ip) is dominated by those 
graphs Go with largest Zq so that, 

Hip -> 00) - p^ Go (18) 

and it is therefore of some interest to identify these maximal graphs. This identification alone is not 
sufficient to solve the models at large p because the number of such graphs and fluctuations around them 
must be included to find the sub-asymptotic behaviour and, in particular, the value of 7 s tr- 

In the following section we identify most of the maximal graphs for models I, II and III. For a given 
value of (3 we might suppose that there is a value p c (P) of p above which the maximal graphs dominate. 
(Of course it may be that the transition to dominance by maximal graphs is seamless and that there is no 
critical value or that p c {P) — 00 for all (3.) In section |] we begin to address the question of the behaviour 

0fp e (/3). 
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3.2 Model I 



In this section we prove that the maximal graphs for model I are tree graphs with tadpoles at the ends of 
the branches (see fig[l] for an example). 

Starting with an arbitrary N vertex </> 3 graph, select a point D and reconnect the links to it as shown 
in fig|2[ The resulting graph still has N vertices, and provided that the link DC lay on a closed loop in the 
original graph, the new graph will still be connected. Repeat this procedure taking care not to disconnect 
the graph. Eventually there are no further links that can be cut and the graph is tree-like (ie there are no 
closed loops except for the tadpoles at the ends of the branches). Before the link DC is cut the partition 
function is 

Z\ = ^ Z(S a ,S b ,S c ) exp(3(S a S d + S b S d + S c S d ) (19) 

S a SbS c Sd 

= Z(S a ,S b ,S c )(l + t 2 (S a S b + S b S c + S c S a )), (20) 

SaS b S c 

where Z(S a , S b , S c ) (> 0) represents the partition function for the remainder of the graph with boundary 
spins S a , S b , S c and C — cosh/?, t = tanh/3. After reconnecting the links, the new partition function is 

Z 2 = ]T Z(S a ,S b ,S c )(l + t)(l+tS a S b ) (21) 

s a s b s c 

so that, 

Z 2 -Z 1= Z(S a ,S b ,S c )t(l + S a S b )(l-tS b S c )>0. (22) 

Thus at each step the partition function increases (equality occurs for — and (3 — oo). The partition 
function takes the same value, 

Z tree = (l + tanh/3)* +1 (23) 

for all tree graphs with TV vertices, so we have proved that the set of tree graphs is maximal for all finite 
non-zero (3. This set of graphs does not magnetize at finite 0. We will need later the number of tree 




Figure 2: Partition functions 
D 





Z(S a ,S b ,S c ) 



Zl 




graphs with N vertices which is given by 



Gtree(N) 



(N-2)\ 



(|JV + l)!(iJV-l)! «- 



(24) 



In considering the contributions of different graphs to Zn(j>) at finite p (which we will do in sections || 
and @) it is useful to examine the ratio, 



1 = "- 



(1 + t) 2 



(1 + S a Sb — t(SbS c + S a S c )) G ^ 



where 



(Q) 



z 2 



J2 Z(S a ,Sb,S c )(l + t)(l + tS a S b )Q. 



(25) 



(26) 



S a S b S c 



For small t, (S a Sb) — 1 or is of order t (depending on whether or not A and B are distinct points), while 
(SbS c ), (S a S c ) ~ 0(t m ) with m > 1; thus Z 2 jZ\ increases with t at small t so that tree-like graphs are 
becoming more important. Assuming that p c (P) is finite, it is decreasing with /3 in this region. On the 
other hand, for large enough t, {S a Sb), (SbS c ), (S a S c ) ~ 1 and Z 2 /Zi decreases towards one as t — ► 1, so 
that tree-like graphs become less important again. The position of any minimum of the curve p c (P) must 
lie between these two regimes. 



3.3 Model II 

For model II, the maximal graphs are ring graphs (fig |3|) and again this is true for any value of (3. The 
proof, which is very similar to that for model I, consists of two parts. Firstly, we show that any graph in 
this model can be converted into a ring graph by a series of steps, where none of the intermediate graphs 
contain tadpoles. Secondly, we show that the partition function increases at each step and thus that ring 
graphs are maximal for all (3. 



G 



Figure 3: Ring graph 




For 4> 3 graphs of genus zero the number of faces, F, is related to the number of edges, E, through 
F = 2 + ^E. Defining fi to be the number of faces with i sides, then F = ^2 fi an d ^ = 5 E so that, 

E^ 1 -^ 2 - ( 27 ) 

i 

Since fi > 0, in order for the equation to be satisfied we must have fi ^ for some i < 6. Thus, any 4> 3 
graph in this model must contain a 2-loop (ie a loop of length two), a triangle, a square or a pentagon. 

Starting with a graph G, containing no tadpoles, replace any subgraphs such as that in figure ^ (which 
we will call dressed propagators) with bare propagators (fig^jb), yielding a new (j> 3 graph, G", to which 
the above theorem applies. Putting the dressed propagators back we recover G and have shown that it 
contains at least one of the following: a dressed 2-loop, or a dressed or bare triangle, square or pentagon. 
The only exception is the ring graph, which upon making the replacement shown in fig ||just gives a circle 
for graph G' . Thus, any graph, except for the ring, contains one of the subgraphs in the above list. 

Dressed propagators containing n 2-loops will be drawn as in fig |^c. We are going to replace the 
subgraphs fig [|a to fig ||a with subgraphs fig ||b to fig ||b respectively. The number of vertices is unchanged 
and the number of 2-loops is increased by these replacements. In appendix [A| we show that the partition 
function increases, for any choice of dressed propagators in the original subgraph (ie for all choices of 
j, k, l,m,n> 0; for the 2-loop case, fig [|a, n, m are not both zero). Thus by choosing the orientation of 
the replacement, we can create a new graph, which is connected, has no tadpoles, has the same number 
of vertices as the original and has a larger partition function. By repeatedly eliminating the subgraphs in 
the above list we will eventually end up with a ring graph, proving that the ring graphs are maximal for 
all /3. The partition function for the ring graph is given by, 

Z„ng - (l+t 2 ) f + (2t 2 ) f (28) 

and this graph does not magnetize for any finite (3. 
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Figure 7: (a) Square (b) Replacement subgraph 




d m c 



Figure 8: (a) Pentagon (b) Replacement subgraph 
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3.4 Model III 



In this case we have not found a procedure which converts any given graph into a maximal graph; essentially 
this is because the form of the maximal graph now depends on (3. We have identified the maximal graphs 
for the limits of small and large f3, but not for intermediate values of /3. 

A high temperature expansion (/3 — > 0) of the partition function for a given graph gives, 

Z G = l + Y j nit\ (29) 
l 

where n\ is the number of closed, but possibly disconnected, non-self-intersecting loops in the graph that 
contain I links. For model III, m = ri2 = because we have eliminated tadpoles and self-energies and the 
first non-zero coefficient is n 3 , the number of triangular loops. The maximum possible value of n 3 is the 
integer part of y which we will denote by [y] . Suppose that TV is divisible by three, then taking any graph 
with points and replacing each point with a triangle (fig |j]a), gives an N vertex graph with n 3 = y . 
If N is not divisible by three, replace all except one or two of the points giving n 3 = [y] ■ To show that 
this really is the maximum possible number, first note that the only case which has two triangles back to 
back is the tetrahedron (fig|^b). Then, ignoring this case, given an N vertex graph we can collapse all of 
its triangular loops to points and get a graph with AT' vertices and no tadpoles or self-energies. Clearly 
N' > Ti3 (since each collapsed triangle yields a vertex) and N' = N — 2n 3 (as collapsing a triangle removes 
two vertices). Thus N — 2?i 3 > n 3 and so n 3 < -j. Hence the maximum value is n 3 = [-3] for N > 4. 



Figure 9: (a) Replacement (b) Tetrahedron 



(a) 




(b) 



Consider graphs, G, for which the total number of points is N — 2 x 3™, where n is some sufficiently 
large integer. In order to maximize Zq in the /3 — > limit, we maximize each coefficient in turn. Having 
maximized n 3 as described above by replacing each point of G' (with |iV points) with a triangle we now 
choose G' in order to maximize the next coefficients. The replacement of points with triangles doubles 
the number of edges bordering each face of G' , which being a model III graph has a smallest loop length 
of three. Any such loops will be doubled to length six; thus n 4 = n 5 = 0. Now, n e = n% + |n 3 (n 3 — 1) 
(where n\ is the number of connected loops of length six), so that we need to maximize n\ next. To do 
this we must maximize n 3 of G' (since loops of length three in G' become those of length six in G) . So we 
take a graph with j-N points and make the replacement of points with triangles to get G' . Carrying on 
in this fashion we end up with a fractal graph (see fig [l0|). If is not a power of three, then the graph 
will not quite be regular, but will still be essentially fractal-like. 

At large (3, graphs which magnetize can be studied in the low temperature expansion for which, 
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Figure 10: Fractal graph 




Zg = Yo 2 ^' 3 y 1 + £ mrX ) ~ % exp N + £ flsa;S ) ' (30) 

where a; = e~ 2/3 and m r is the number of domain boundaries that cross r links. For graphs with no 
tadpoles or self-energies any domain boundary must cross at least three links so that the sum starts at 
r = 3 and the sum for the free energy also starts at s = 3. However, this does not apply to graphs that 
do not magnetize; the ladder graph (figO) has m 3 oc N, but 7714 = ~ — l) so we might expect it to 
exponentiate to give a series starting at s = 2 with 02 = \ . It is straightforward to check that this is the 
case from the ladder free energy, which is given by, 



\ log (\ (l + f 2 + v/(l- 42)2+4*4)) . (31) 



gladder 

Because these graphs are one-particle irreducible it is not possible to get 7712 oc iV 2 and a\ is always 
zero. We conclude that at large (3 the ladder graphs, whose free energy starts at 0(x 2 ), will dominate 
magnetizable graphs, whose free energy starts at 0(x 3 ). Furthermore it is easy to check that gadder is 
bigger than that of the fractal graph for large (3. 

What happens at intermediate (3 is not clear but is probably quite complicated, involving graphs which 
in some sense interpolate between the fractal and the ladder graphs; whether the transition from fractal 
to ladder is continuous or discontinuous we cannot say, but the former seems more likely on grounds of 
universality with models I and II - but maybe this system is not universal. Neither the fractal nor the 
ladder magnetize so probably the intermediate graphs do not magnetize either and as p — > 00 the model 
only magnetizes at f3 — > 00. 

It is interesting to note that for a Gaussian model on a random triangulation embedded in D dimen- 
sions |27|, |2^, |2^, ||(], |3jJ the maximal graphs (for D > 0) are those with the minimal number of spanning 
trees and that the corresponding </> 3 graphs are tree graphs, ring graphs and ladders for models I to III 
respectively. Thus the maximal Ising graphs at large (3 are the same as the maximal graphs for the Gaus- 
sian model, but they apparently differ at small (3 in the case of model III. Reference |27| lists the maximal 
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Figure 11: Ladder graph 




graphs, but gives a more symmetric version of the ladder graph; presumably that is because fig |llj is 
excluded from the relevant model even though it has no tadpoles or self-energies. 



4 Truncated models 



4.1 Model I 

To get some indication of how large p must be before the maximal graphs become dominant we now 
consider a truncated model in the weak coupling regime. For small /3, 



1 

JV 



1 



rut + I n 2 - -iii 2 )t 2 + [ n 3 - mn 2 + -ni J I t 3 + 



1 



3 \ .3 



so we consider a model for which fie is truncated to 



t n i . 



The grand canonical partition function for this model is then, 

oo oo 

Zfapt) = J2 e '" N E ^ (1) (^^i)e" ipt , 



(32) 



(33) 



(34) 



ni=0 



where G W (N, n\) is the number of graphs with N vertices and n\ tadpoles (loops of length one) 
(0 < n\ < + 1). We show in appendix |B.2| that G^{N, n\) satisfies the recurrence relation, 



0« (JV + 2, m + 1) = ( 3N 2 ? x ) gW (JV, nt) + 2Q& (JV, m + 1) 
\ ni + 1 / 

and the total number of graphs is known so 



(1) (iV,nx)=aW(JV), 



(35) 



(36) 



rii=0 
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where Q^(N) is given in (7). Using the recurrence relation (35) and putting y — e pt , x 
that Z satisfies the differential equation, 

^(l + 2.^-l)) = 3x^ + ^, 



e M wc find 



(37) 



which has the solution, 



1 1 i 00 

Z = (hi l) + - (y 1) - - log/* + G {1 \N) x N h- 



(38) 



where h = (1 — 4a; 2 (y — 1)). The asymptotic behaviour at large iV is thus given by 



£~^5«(JV)A-t -J]iV-5exp-Ar('/i+^log^-ilog(12V3)") , 

AT iV ' 



(39) 



so the thermodynamic free energy /i c (pt) obeys the cubic equation, 



/i c + log (1 - Ae- 2 ^(y - 1)) - X - log(12V3) = 0. 



(40) 



The solution is 



/, ' ;/ " : log2 ' \ log(y 11 ' - — " I W^W 



E 

<T=±1 



18 



CO 



+ 



54 



(y-i) 2 



(41) 



with w = for y < 3 and u> = 1 for y > 3 (the solution and its derivatives are actually continuous 
across y = 3). The solution is plotted in fig 12a where it is compared to fi c ' = \pt + log 2, which is what 
we would expect if trees were totally dominant. We see that n c {pt) approaches fi c ' quite rapidly as pt 
increases, but it is necessary to have p — > oo for the trees to dominate completely. It is also interesting 
to examine the ratio r = (ni)/N, which is shown in fig 12b; for pt > 6 the ratio is very close to one half 
and the fact that the gradient changes very quickly near this point leads us to suspect that the full model 
might develop a discontinuity in one of its derivatives at some finite value of p (ie p c (/?)). A formula for r 



46^ 



can be derived from that for Z; putting b = 2e pt , A = ( 3fc _ 2 )(b + 2) > 

1 



2-6 



1 + . 



aj L(i- VT~^a) 3 + u? (i + VT^xj ' 



(42) 



where w is a cube root of unity and we need to use w = lfor0<fo<| and w = e^ 1 for | < b < 2 (the 
solution and derivatives are continuous across b = |). 

As is shown by (39) the string susceptibility "f str = —\ for all finite p in this truncated model; only 
for p = oo does 7 str change to |, but, again, in the full model this change may occur at some finite p c {0). 
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Figure 12: Truncated model (model I) 
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4.2 Model III 

This calculation can be repeated for model III, which has, 

1 



fJ-G = T7(M 3 + n 4 t 4 + ri 5 t 5 + (n 6 - -n 2 3 )t 6 + ■■■). 



N 



We truncate /ig to 
where < n 3 < [y] . The grand canonical partition function is defined as 



Mg- N t 



(43) 
(44) 

(45) 



Z(fji,pt 3 ) = E e ~" N E 5 (3) W«3)e n3pt , 

N=6 7i. 3 =0 

where g^(N, 713) is the number of graphs in model III with N vertices and 713 loops of length three. This 

(46) 



obeys a recurrence relation, derived in appendix JB.l 

G {3) (N + 2,n 3 + l) 



N - 3n 3 



n 3 + 1 



g^(N,n 3 ) + 3g^(N,n 3 + l) 



(3), 



which holds for iV > 6. The recurrence relation (46) gives, putting y = e pt and a; = e 

92 - „ , , . \ o dZ 1 ft 
_(l + 3^ (y -l))=x3_ + r V 



The solution is 



z( M ,^ 3 ) = -i (hx - \x> + \A + £ g^(N)h N , 



where 5 (3) (^) 

is given by (9) and h = x + x 3 (y — 1). Thus for large N, 

Z-Y. S (3) (N)h N ~ E N-i expiV (\ log (^) + log (a + X 3 (y 1)) 



A' 



A' 



(47) 
(48) 

(49) 



Solving for /j c (pt 3 ) as in the previous case gives 



1 



fi c {pt 6 ) = --log3 + -log2+-log(y-l)-log 



2 3 
We can also calculate (n 3 )/N, 

(«3> 1 



1 + °V 1 + (I) iPT 



cr=±l 



(50) 



TV a 



x _ 1 / 3-a 



2 V 12 - ca 



E 1+^/1-4 



<X=±1 



3 — a 
12- ca 



(51) 
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where a = 3(1 - e~ pt ) and c = 295/256. A gain, this saturates as p is increased (see fig 13b). 

Looking at the graph (fig 13a) of [i c (pt 3 ) — ^pt 3 — ilog(^) (where the last term takes account of 
the exponential number of graphs with = |) we see that for pt 3 > 8 the behaviour is essentially linear, 
because after this point the integral is being dominated by a single type of graph (ie those graphs with 
n 3 rts jjN). However, it should be noted that in this truncated model we again need p = oo in order to 
actually reach n^/N = -| and the set of such graphs is quite large (not just fractal-like graphs). This is 
due to the fact that \iq has been truncated to the first term (n.3), but higher order terms are needed to 
show that fractals are maximal for j3 — > 0. Again, ^jstr — — o 

for all finite values of p. 



Figure 13: Truncated model (model III) 




5 Magnetization transition 
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5.1 Derivation of a bound on (3 C 

In this section we derive a bound on the critical value of the coupling constant, /3 C , for the case of a single 
Ising spin (p = 1) on a fixed graph G. This implies a bound on the critical coupling for magnetization 
related phase transitions in models I, II and III. The high temperature expansion of Zq (29) gives, 

Z G = 1 + Y,nit l =explf^ ai A. (52) 
i \i=i J 

Now consider, 

Z G '=exp^n?t^, (53) 

where n\ is the number of connected, non-backtracking closed loops of length I. Expanding the exponential 
yields all closed non-intersecting loops, but in addition it gives loops in which some links are used more 
than once so that Zq < Zq> and therefore 



N 
i=i 

However the number of non-backtracking closed loops of length I originating at a given point is less than 
2 l so 

OO 

i=i 

Thus the high temperature series for hg must converge if t < ^ ; it follows that any phase transition on G 
must occur at a t c > \ (ic f3 c > 0.549). 

A slight modification of this argument excludes the possibility of permanently magnetized graphs. Fix 
one spin S + to be +1, so the magnetization Mq(0) is given by, 

+ {S} i (a, 6) 

where the runs over au possible locations of the fixed spin (this is necessary because unlike the regular 
+ 

lattice model any two given spins may be almost disconnected from one another). The only contribution 
to the numerator comes from paths connecting Si to S + , so 

E 2K "'EE*(«. +)«' 

MoOT " s* + s n ' <57) 

{S} <o,6> 

where d/(i,+) is the number of paths of length I from ^ to S+ (which can be disconnected, but are 
non-self- intersecting and non-backtracking). However, 

Y;di{i,+)t l =^wj(i,+)t ( (l + feif + & 2 f 2 + •••), (58) 
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where wi(i,+) is the number of connected paths of length I from Si to S+ and the bi series gives the 
contributions from the closed loops, which do not intercept the path. The denominator, however, contains 
contributions from all closed loops and hence is larger than 2 JV_1 J^i bit 1 so 

m g^<^EEE^( 8 .+)*' ( 59 ) 
+ t / 

However, J}» Wi(i, +) is just the number of connected paths of length I from S+ and we know that this is 
less than 2 l . Thus, 

Mc^X^E^l^, (60) 

1 

for t < \ and hence Mq(JS) — > as N — > oo. If more than one spin is fixed to be +1, then the right hand 
side of this equation is multiplied by the number of such spins, so that to avoid Mq{(3) — > a number of 
spins proportional to N must be fixed. That is, we need to fix a thermodynamically significant number of 
spins in order to get the system to magnetize and thus there can be no spontaneously magnetized state 
for t < \. 

This result was derived for a single graph, but obviously still applies when we perform a summation 
over graphs; no phenomena associated with magnetization can occur at t < | . 

For the </> 3 graphs the high temperature (f3 — ► 0) expansion converges for tanh/3 < |, this implies 
that on the dual triangulation the low temperature (/3 — + oo) expansion converges for e~ 2/3 < i (ie for 
(3 > | log 2 f» 0.347). Thus any critical value of (3 must satisfy j3 c < \ log 2, on the triangulated surface. In 
contrast with (f> 3 graphs, there exist both graphs which never magnetize and those which are permanently 
magnetized (see fig 14). The argument showing that there are no permanently magnetized (f> 3 graphs, can 
not be used for the dual triangulation because the relation ^ i Wi(i, +) < 2 l no longer holds. This is due 
to the fact that the vertices can have coordination numbers greater than three (in fact the permanently 
magnetized graph shown has points with coordination numbers of order N). 



Figure 14: Triangulations: (a) never magnetizes; (b) permanently magnetized 



(a) (b) 
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5.2 Mechanisms of magnetization 

For this model the magnetization can be written as 



E te m g{z g ) p 

M(p,/3)= — = r -— , (61) 



G 



where Mg{(3) is the magnetization for a given TV-vertex graph G, and we are taking the limit N —* oo. 
In the limit p — > we effectively have a quenched magnetization, 

M (f3) = M(0, 0) = -Lj £ -^- m gW), (62) 

where Q(N) is the number of graphs with N vertices in whichever model is being looking at. This case 
has recently been investigated numerically [f32[ . The critical value for the magnetization of this model will 
be denoted by (3 C and the critical coupling constant for a given graph G by /3q. One can show that the 
behaviour of Mo(/3) at /3 C only depends on those graphs which magnetize near (3 C . We will illustrate this by 
proving the result for the case in which graphs G undergo a first order phase transition; in this case only 
those graphs for which j3 c — e < (3q < /3 C , (where e is an arbitrary small positive number) contribute to 
Mo(flc). In actual fact the phase transitions of individual graphs are second order, however the extension 
of the result to this case is relatively straightforward and is omitted. Now, 



M„(/3) 



g(N) 



Y M G (fi) + Y M g(/5) 

G:l3^<fj-e G:0^>l3-e 



(63) 



where the notation G : 0q < (3 — e means that we are summing over graphs G for which the inequality is 
satisfied and we have absorbed the symmetry factors into the summations. Since the system magnetizes 
at p C) 

M (p c -e) = -^-r Yl M o(Pc-e) ^ 0. (64) 
For a first order transition in which each graph jumps to Mq at its critical point, 

m^. MoV '- e) - mJ^J? - mm(M ° a) sm S 1 ■ («) 



While, 



which tends to zero by (64). Hence (63) gives 



M (/3 C ) = ^- E M o(Pc) = ^ E M ^o) , (67) 

y > G:p G >(3 c -e VK > G:/3 c - £ </3J,</3 c 
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so that only graphs which magnetize near (3 C contribute to Mq((3 c ). As it stands the proof cannot be used 
for second order transitions (where M<2, = 0), but careful consideration shows that in this case only graphs 
for which f3 c ~ e < /3q < (3 C + S (where e, 6 are small positive numbers) can contribute to Mq((3 c + S). 
Thus the critical exponents can only depend on the behaviour of graphs which are magnetizing near the 
critical point f3 c and the magnetization is as a result of spin-ordering of these graphs. 

For non-zero p, a different geometric type of transition is possible and this seems to be what occurs 
for large p. In this case the magnetization is caused by the changing of the relative weights between 
different graphs. This can best be understood by looking at a simple model in which there are only two 
types of graph; suppose that there are n\ unmagnetized graphs (ni ~ exp f\N) with partition function 
Z\ ~ exp/iiA^ and n 2 magnetizable graphs (r%2 ~ exp f%N) with magnetization m{[3) and partition 
function Z 2 ~ exp [i 2 N . Then assuming that the magnetizable graphs arc the more numerous {f 2 > fx) 
and have smaller partition functions (fi 2 < fJ-i ) , which is certainly the case if we are looking at transitions 
to tree graphs in model I, the partition function is 

Z N {p) = n x Z x v + n 2 Z 2 p ~ e (/i+PMi)JV + e (/ 2 +PM2)w ( 68 ) 

N— .00 

and the magnetization, 

M{p,(i) ~ — Vt m(S) e ^ 2+p ^ N . (69) 
Zn(P) 

In the thermodynamic limit the magnetization M(p, (3) is zero if fx+pfJ>i > /2+PM2> that is, if p{^x — ^2) > 
(/2 — f\). As (3 is varied the difference \i\ — fi 2 changes and at the point for which the inequality no longer 
holds (/3 C ) the magnetization jumps from zero to m{(3 c ). Thus there is a first order transition. 

In the large p limit for which /3 C — 00 it is possible in model I to have magnetized graphs with r = m/N 
arbitrarily close to i, so that even if the transition is generically first order the discontinuity, Ar, could 
be zero in this limit (ie perhaps Ar — > as p — > 00). For p — > 00 there is such a transition, for model I, 
between unmagnetized tree graphs and magnetized non-tree-like graphs, with no discontinuity fHj, 16 . 



6 Conclusion 

In figure |l5| we have drawn possible forms of the phase diagram for model I. The magnetized region 
is labelled M, the tree-like region T and the remaining unmagnetized non-tree-like region U. We would 
expect similar diagrams for models II and III, but in the latter case it is not clear what is happening in 
the T region, as we have not identified the maximal graphs for intermediate values of /3. 

For models I and II, we know that the partition functions of graphs increase as we replace closed 
loops with disconnected structures. Since the graphs are becoming less well-connected one might expect 
that they tend to magnetize at higher values of (3 and have lower magnetizations at a given value of j3. 
As p is increased in this model, the graphs with the larger partition functions make a relatively larger 
contribution to Zn(j>) and so we expect that (3 c (p) will increase. There is some evidence for this from 
numerical simulations [jl9| pl[ ^3). We might also expect the magnetization, for a fixed (3, to decrease 
as p is increased and again there is some evidence for this from Monte Carlo simulations |^| . It should 
of course be noted that these simulations were for model III, for which we do not know the maximal 
graphs. However, since all the maximal graphs that we have identified are not well-connected and do not 
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Figure 15: Possible phase diagrams 
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magnetize, it is quite likely that this also applies for model III at intermediate values of /3; it also follows 
that in none of the models do we expect (3 c (p—KX)) to be finite. 



In the papers of Wexler Jig , 16 the p — - > oo limit of model I is studied in the low temperature (large 
P) expansion. In this limit there is a critical coupling e _2/3c ~ 1/p below which the dominant graphs are 
tree-like being made up of large numbers of almost disconnected baby universes within which the spins are 
all aligned. This phase has 7 str = i and our results indicate that this behaviour must persist beyond the 
large j3 expansion, in fact for all f3 > 0. Above j3 c the dominant graphs contain few baby universes each of 
whose volume (number of vertices) is very large and 7 str = — \\ thus it seems that the properties of the 
M region do not depend strongly on p. At the critical point, where both the number of baby universes 
and their volume diverge, it is found that "f s t r = This result is also found in the paper of Ambj0rn et 
al fT^I who considered a restricted model in which only phase interfaces of minimum length are allowed; 
they obtain the result that the exponent at the critical point 7* is related to that in the large (3 phase, 
7, by 7* = 7/(7 — 1) so long as 7 < 0. This phase transition at p — > 00 is geometrical, as discussed in 



section 5.2, rather than one of spin alignment; it occurs because the partition functions of magnetized 
graphs are catching up with the others. It seems highly likely that this phenomenon will persist to finite 
p although a calculation of the 1 /p corrections (or equivalently the contribution of domain boundaries of 
length greater than two in the model of [ fl7| ) is needed to rule out the phase diagram of [l5[b analytically. 

The extent of the tree-region T is not yet fully understood. We know that at (3 = 0, where the Ising 
partition function is one, the model is independent of p and has jstr = —\\ in addition, the absence of 
extra transitions in models with c < 1 suggests that the line AC separating the U and T regions does not 
go below p = 2. We cannot rule out the possibility that p c (P) — 00 as shown in figure |l5|a. However, the 
truncated models do display a region of relatively rapid changeover to the T region and it seems likely 
that the full models have a phase transition (with p c t ps 6, p c t 3 w 8 for models I and III respectively at 
small enough t). 

For model I, x = | — r (where r — ni/N is the ratio of loops of length one to vertices) serves as an 
order parameter. In the region T, x — and approaching the line AC from below at fixed /3, we would 
expect x to decrease continuously to zero, that is a second (or higher) order phase transition. This is 
because, just below the critical line, the largest contributions to Zn(p) come from diagrams which are 
tree-like except for a number of closed loops. As p is increased the number of closed loops (excluding 
loops of length one) decreases to zero in a continuous fashion (see for example fig 12b, which shows how 
(r) approaches a half as p is varied, for the truncated model). 

The transition to T from M across CB is rather different. Then we know that just below the critical line 
the partition function is dominated by magnetized graphs (which cannot be tree-like or almost tree-like, 
as such graphs do not magnetize for t 7^ 1) and such graphs have r 7^ i. Thus we might expect there to 
be a first order transition with r jumping discontinuously. It is tempting to hypothesise that for the line 
DC between the M and U regions there is a continuous spin-ordering transition whose characteristics are 
determined by the graphs which are actually magnetizing at the critical point /3 c (p), as happens for the 
quenched case (p — > 0). On the other hand the line CB between the T and M regions is a geometrical 
phase transition in which the transition is caused by the changes in geometry due to the changing relative 
weights (ie Zg p ) of the magnetized and tree-like graphs. 

For model III, there is direct numerical evidence from computer simulations. In this case it is not as 
clear what we should use as an order parameter for the U T transition, but n^/N (the ratio of closed 
loops of length five to vertices) seems a plausible candidate. For both the fractal and ladder graphs this 
quantity is zero and if, as seems possible, the intermediate maximal graphs consist of mixtures of fractals 
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arranged in a ladder-like fashion, then it would be zero in the whole T region. Unfortunately, none of 
the numerical studies looked at this quantity in detail. Reference Q measure the quantity for p = 16 
(in fig 12 of that paper), but only at the critical point /3 C , where it is very close to its pure gravity value. 
However, references M, do give graphs of n^/N for various values of (3. This quantity peaks at 

some value of (3 below (3 C and drops to the pure gravity value as (3 — > or (3 —* (3 C . The height of the peak 
grows with p in an apparently linear fashion, but we know that it is bounded by the relation n 3 /N < i. 
It would be interesting to see how this quantity saturates as p is increased, for (3 near /3 C , and to compare 
it with the results we have from the truncated model in the limit of small (3. 

It is interesting to examine the average maximal radius, r max , of the ensemble. For a given graph r max 
is defined in the following fashion. Take an arbitrary point on the graph (called the "centre"), mark all of 
its neighbours as being at a distance one from the centre, mark all the unmarked neighbours of these points 
as being at distance two and so on. The maximal extension is then the distance of the furthest point from 
the centre and averaging over different centres gives r max on the graph. This procedure can be carried 
out on either the <fi 3 graph or the dual triangulation, giving two different definitions, which we will denote 
r max an d r max respectively. In figure 7 of pp| r^ lax is plotted for various different models; the curve has a 
minimum for values of (3 just below (3 C and its depth increases with p up to the maximum value considered, 
p = 16. Both fractals and ladders have very small values of r^ aa . and so as the T phase is approached one 
would expect r^ aa . to fall. Let p* be the value of p at the point C in fig |l5| c and now suppose that p is 
held fixed at p < p* and we approach the magnetization transition from the U region; then as (3 increases 
we also get closer to the T region and thus might expect to see a decrease in r^ ax , which disappears 
when the system magnetizes. As p | P* this effect would become more pronounced until at p just above 
p* we expect to see a sharp transition into a branched polymer-like phase which immediately disappears 
as (3 increases and the system magnetizes; r^ lax would decrease sharply as the critical region around C is 
approached and then rebound when the system magnetizes. This behaviour seems very like that observed 
in pC| ] which suggests to us that even the largest value of p used in the simulation is probably no greater 
than p* . If p is significantly greater than p* we would expect the minimum value of r^^, when the 
ensemble is most branched polymer-like, to occur somewhere in the middle of the T region rather than 
at the magnetization transition; so as p increases beyond p* the minimum of r^ aa , will move away from 
the magnetization transition. It is difficult to be entirely certain in interpreting the behaviour of r^ aa . 
without knowing how the minimum value scales with TV; it would be useful to have simulation data on 
this and also at much larger p to be certain of our interpretation. 

Although we have worked entirely with Ising models many of our considerations probably extend to 
other spin systems which have second order phase transitions on regular lattices; in particular the maximal 
graphs and the general behaviour of the truncated models are likely to be the same. 

We would like to thank Thordur Jonsson for giving us advanced details of [|l7j and to acknowledge the 
support of the SERC under grant GR/J21354 and through the research studentship GR/H01243. 



A Proof that ring graphs are maximal 

In order to complete the proof that rings are maximal for model II, we need to show that each of the 



subgraphs listed in section 3.3 can be eliminated in a way that increases the partition function. 
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Consider the case for the 2-loop first. A propagator dressed with n 2-loops and joining spins S a to Sb 
(fig [|c), has a contribution of 

G (n) {S a , S b ) = 2 2n C 3n+1 ((1 + t 2 ) n + S a S b t 2n+1 2 n ) = C'(l + t'S a Sb), (70) 

where C = cosh (3, t = tanh (3 and 

C' = 2 2n C 3n+1 (l + t 2 ) n (71) 
/ 2t 2 \ n 

'' °° (it?) ' s Y " L (72) 

Thus dressing a propagator effectively renormalizes the coupling constant for that link. The partition 
function for the whole of the original graph is given by 

Z(S a ,S b ) Zi(S a) S b ), (73) 

s a s b 

where Z{S a , Sb) is the partition function of the remainder of the graph (Z(S a , Sb) > 0) and Z\ is that for 
the subgraph being replaced (fig ||a). The partition function for the graph after the replacement has Z\ 
replaced by Z% (fig |5|c). The powers of two and the factors of C and C" are the same for both graphs and 
cancel with terms in Zq , so that we can drop these factors which gives, 

Z 1 = l + Y n+m t 2 + S a S b t 3 (Y n + Y m ), (74) 

Z 2 = 1 + t 2 + S a S b (2t 3 Y n+m ) . (75) 

Hence, 

AZ = t 2 (l-Y n+m + S a S b t{2Y n+m -Y n -Y m )) (76) 

= t 2 {{l-Y n ){\- S a S b tY m ) + Y n {l~Y m ){\~ S a S b t)) >0. (77) 

Noting that < Y, t < 1 and n, m > 0, we can see that AZ is not negative and that AZ = occurs for 
t = and t = 1, as we would expect. The partition function of the graph is increased by the replacement 
in fig H 

The triangular case (fig |6|) can be proven in the same way as the 2-loop case. However, for the square 
and pentagon, the expression for AZ is so complicated that it is extremely difficult to regroup the terms 
and write AZ in a form which is manifestly non-negative. An inductive method can be used to prove each 
of the cases and we will indicate below how one would use it to prove the case of the square (fig [?]) . For 
this case we need to prove that AZkimn{S a , Sb, S c , Sd) where 

AZ k i mn t- 2 = 2+t (S a S b Y k (2Y" - 1) + S c S d Y m (2Y l - l) - S b S c Y l - S d S a Y n ) + 
t 2 (1 - Y k+l+m+n - S a S c (Y k+l + Y n+m ) - S b S d (Y m+l + Y n+k )) + 
t 3 (S a S b Y n (2Y k - Y l+m ) + S c S d Y l (2Y m - Y n+k ) - S b S c Y k+n+m - S d S a Y m+l+k ) + 

t 4 S a S b S c S d (4Y k+l+m+n - Y k+m - Y n+l ) , (78) 
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is positive for all < t < 1, and all values of the spins, S a , Sb, S c , Sd- It is easy to prove that AZoooo > 
for < t < 1 and any values of the spins S a , Sb, S c , Sd- We define a difference operator Dk by 

DkAZklmn = AZ( k +l)lrnn — YAZklmn (79) 

and similarly for Di, D m , D n . Now, AZkimn has the form, 

AZumn =A + Y h B, (80) 
where A, B are independent of k, but in general depend on all the other variables. Hence, 

D k AZ Hmn = (1 - Y)A, (81) 
which is independent of k; note that (1 — Y) is a positive factor. Now, 

AZ( k +l)lmn = YAZklmn + DkAZklmn- (82) 

If we can show that Dk AZkimn > for any Z, m, n and similarly for D m and _D„ acting on AZ , then, 
since AZoooo > 0, the result AZkimn > for all k, I, m, n > follows by induction. 

In order to prove that DkAZklmn > we use the fact that the difference operators commute so that, 

DkAZ u (m+l)n Y(DkAZklmn) + Dm(DkAZklmn)- (83) 

Now we need to prove that DkAZkoao > and that D m DkAZkimn > (and similarly for all pairs 
of difference operators). Proceeding in this way the problem is reduced to showing that the following 
quantities are positive: AZoooo, DkAZkooo (similarly for Di etc), DkDiAZkioo (similarly for all possible 
pairs), DkDiD m AZkimO (all possible triples) and DkDiD m D n AZkimn- Each of these expressions depends 
only on t and the spins, and is thus relatively easy to prove (especially if a computer is used to check each 
case). Thus finally we have AZkimn > for < t < 1. The pentagonal case can be proved in the same 
fashion. 



B Derivation of recurrence relations 
B.l Model III 

The recurrence relation (46) for two-particle irreducible planar graphs is derived by considering the effect 
of replacing a point P, on a graph G, by a triangle. Such a replacement increases the number of vertices in 
the graph by two and increases the number of triangular loops (713) by either one (figj^a) or zero (fig |l6|) 
depending on whether or not the point P of graph G was lying on a triangle. Note that since self-energy 
terms have been eliminated the only possible diagram for which P could lie on two triangles is given by 
fig||b, which has N — 4 vertices; we shall ignore this case and derive a recurrence relation valid for N > 6. 

Graphs with N + 2 vertices and ^3 + 1 triangles can be made from those with N vertices and 
triangles by replacing one of N — Sn^ points (fig ||a), or from graphs with N vertices and + 1 triangles 
by replacing one of 3(713 + 1) points (fig |l6| ). We will refer to the points which yield the correct number 
of triangles when replaced (ie 77,3 + 1) as being 'allowed' points. Clearly all such graphs can be made by 
this method. 
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We need to determine the effect of such a replacement on the symmetry factor of a graph. Suppose 
that graph G has symmetry group Q of order s. Since G is planar we can think of it as being a polyhedron 
and the members of Q as being rotations. By making rotations we can interchange a point P with any of a 
set R (of size /i) of equivalent points. Now, rotations that leave P in the same position give a subgroup Ji 
of Q which has order s/fi (which will equal either one or three, since there are three lines emerging from 
each vertex). Replacing P by a triangle T, we group the triangles of the new graph G' into equivalence 
classes. If T is a member of a class of size / 2 then the order of the symmetry group of G' is s' = sf 2 / fi- 

If we consider the equivalence classes of all the allowed points in all graphs with N vertices, then there 
is a one-to-one correspondence between these classes and the classes of triangles in all the graphs with 
N+2 vertices. A replacement of a point P by a triangle T defines a mapping between classes of points 
and of triangles. Let us associate a weight fi/s = /a/s' for this mapping between classes. That is, we are 
associating a weight of 1/s for each allowed point in a graph that is being mapped from, or equivalently 
a weight of 1/s' for each triangle in a graph that is being mapped to. The total weight for allowed points 
is given by, 

E E 7= E (N-3n 3 )±+ £ 3(n» + l)i (84) 

qraph i allowed graph i graph j J 

points <N,n 3 ) (N,n 3 + 1) 

= (N- 3n 3 ) £ (3) (N, n 3 ) + 3(n 3 + 1) G^(N, n 3 + 1). 

However, this is equal to the weight for the graphs being mapped to, due to the one-to-one correspondence, 
which is 

E E 7= E (n 3 + l)^(n 3 + l)g^(N + 2,n 3 + l). (85) 

graph i trianqles 1 qraph i 1 

(N + 2,„3 + l) 

Equating these gives the claimed result. 

Figure 16: Replacement does not change n 3 




B.2 Model I 

The recurrence relation (35) for model I can be derived in a similar fashion, by considering the addition 
of a tadpole to a link (see fig |l7|a) . There is the added complication that the tadpole can be orientated 
in one of two different directions, (we distinguish between these configurations because they are counted 
separately in the matrix model). We will consider the creation of graphs with N + 2 vertices and rii + 1 
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tadpoles, from N point graphs with either m tadpoles (see fig pT^i) or ni + 1 tadpoles (fig p~7|b) . For the 
N + 2 vertex graphs we will group the tadpoles into equivalence classes. However, there are no rotations 
that leave a given tadpole fixed (since this would correspond to rotating a polyhedron, whilst keeping an 
edge and an adjacent face fixed, which is impossible). Thus, rotations just interchange the tadpoles within 
each equivalence class in a cyclic fashion and the order of the symmetry group, s' , is just equal to the size 
of the equivalence classes (all of them being the same size). (Incidentally, this means that the order of the 
symmetry group of a graph divides the number of tadpoles. Hence, if the number of tadpoles is prime, p, 
then the graph has a symmetry factor of either 1 or 1/p.) 



Figure 17: (a) Replacement (b) Replacement on tadpole (c) Orientated version 



(a) 



(b) 



(c) 
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For the N point graphs, we decompose the graphs into sets of allowed links (that is, links which when 
we add a tadpole to them give us a total of n\ + 1 tadpoles). Each N point graph has ™ links and 3N 
orientated links (ie we count each link twice; once with an arrow going in one direction and once with the 
arrow reversed) . The orientated links are collected into classes which contain links that can be rotated into 
each other. A link may or may not be in the same class as its reversed version. It is impossible to make a 
rotation that leaves a link in the same position and with the same orientation. As a result, we again have 
that the equivalence classes are all of the same size, which is equal to the order of the symmetry group, 
s. Now, the addition to an orientated link (as in fig |l7jc) of a tadpole on the right-hand side, defines a 
mapping between equivalence classes of allowed links in N point graphs and of tadpoles in N + 2 vertex 
graphs. This mapping is again one-to-one. The sizes of the equivalence classes equal the orders of the 
symmetry groups, so that the number of classes for a graph equals the number of allowed links divided 
by the order of the group. Equating the numbers of equivalence classes, we get, 
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(86) 



2 - (1) (tf,m) + 2(7H + l)g<U(N, ni + 1) = (m + 1) G (1) (N + 2,7n + 1), (87) 
which gives the claimed result. 
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